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Summary 


A  spectral  numerical  method  is  used  to  study  the  two-dimensional 
unsteady  flow  over  a  flat  plate  in  the  presence  of  a  plane  sound  wave 
propagating  parallel  to  the  flow.  For  s  =  Wx/U^  >  0(1),  no  obser¬ 
vable  interaction  or  energy  interchange  between  the  sound  wave  and 
the  Tollmien-Schlichting  wave  is  present.  In  the  region  s  <  0(1), 
near  the  leading  edge,  an  interaction  occurs;  the  magnitude  of  this 
interaction  is  computed.  A  plane  sound  wave  with  dimensionless  fre¬ 
quency  =  56  x  10~k  generates  a  Tollmien-Schlichting  wave  of 

the  same  frequency  and  an  amplitude  at  the  first  neutral  point  10~4 
times  the  sound-wave  amplitude.  Computations  at  a  second,  lower 
frequency  result  in  a  lower  neutral  point  amplitude.  The  response  of 
the  boundary  layer  to  large  unsteady  disturbances  is  discussed  and 
compared  with  published  behavior  based  on  the  boundary  layer  equa¬ 
tions.  Solutions  of  the  Orr-Sommerfeld  equation  are  examined  and 
found  to  have  a  different  character  for  s  £  0(1)  consistent  with  the 
more  complete  equations. 

I.  Introduction  r  ‘ 

The  present  work  numerically  investigates  the  flow  over  a  flat  plate 
in  the  presence  of  a  plane  sound  wave  propagating  parallel  to  the  free- 
stream  flow.  It  is  postulated  that  the  sound  wave  generates  Tollmien- 
Schlichting  waves  which,  in  turn,  may  lead  to  transition.  Although 
the  investigation  considers  only  small  amplitude  sound  waves,  the 
numerical  methods  used  are  not  restricted  to  this  case. 

An  experimental  study  of  the  generation  of  Tollmien-Schlichting  waves 
by  sound  waves  has  been  made  by  Shapiro  [  l].  His  data  indicate  that 
a  Tollmien-Schlichting  wave  is  generated  with  the  same  frequency  as, 
and  with  a  much  smaller  wave  length  than,  the  sound  wave.  Thomas 
and  Lekoudis  [2]  compare  Shapiro's  data  with  a  simple  model  in  which 
the  flow  field  is  a  sum  of  a  Tollmien-Schlichting  wave  and  a  sound 
wave.  This  model  is  in  fair  agreement  with  Shapiro  [1];  thus,  they 
conclude  that  there  is  no  interaction  between  the  waves  in  the  range  of 
the  data,  and  that  the  interaction  which  produced  the  Tollmien- 
Schlichting  wave  occurs  near  the  leading  edge. 

The  work  on  unsteady  bounda ry -  laye r  theory  is  relevant  to  the  present 
work,  even  though  the  bounda  ry  -  laye  r  equations  cannot  predict  the 
presence  of  Tollmien-Schlichting  waves.  Early  work  on  unsteady 
boundary- layer  theory  was  done  by  Moore  [3]  and  Ostrach  [4]. 
Illingworth  [5]  has  addressed  the  problem  of  the  interaction  of  a  weak 


sound  wave  with  a  compre s sible „  heat-conducting  boundary  layer. 
He  generates  two  series  solutions  for  the  perturbation  to  the 
Blasius  flow.  One  expansion  is  valid  for  small  dimensionless  fre¬ 
quency,  s  =o»x/U  and  the  other  for  large  s.  Lam  and  Rott  [6] 
consider  a  general  incompressible  boundary  layer  with  a  small 
unsteady  disturbance  at  the  edge.  They  give  special  attention  to  the 
joining  of  the  small  s  and  large  s  solutions.  Ackerberg  and 
Phillips  [7]  consider  in  compressible  flat-plate  flow  with  small  fluc¬ 
tuations  in  the  free-stream  velocity.  Both  investigations  conclude 
that  the  downstream  (large  s)  solution  is  approached  asymptotically 
through  damped  eigenfunctions.  Ackerberg  and  Phillips  present 
numerical  results  to  support  this  conclusion. 

II.  Physical  and  Mathematical  Problem  Description 

The  problem  of  interest  is  the  interaction  of  a  plane  sound  wa  ve 
with  a  Blasius  boundary  layer  in  incompressible  flow.  Thus,  the 
sound  speed  and  wave  length  approach  infinity  and  the  wave  number 
goes  to  zero.  Following  Illingworth  [5],  the  dimensionless  free- 
stream  velocity  is  given  by 


u  =  1  +  e  coscjt  ( 1 ) 

3 

where  e  is  typically  10  in  the  results  reported  herein;  cj  and  t  are, 
respectively,  the  (temporal)  frequency  of  the  sound  wave  and  time. 

Going  to  the  incompressible  limit  allows  a  previously  developed  com¬ 
puter  code  [8]  to  be  applied  with  modifications  only  to  the  boundary 
conditions.  Since  a  detailed  description  of  the  formulation,  the 
assumptions,  and  the  numerics  has  been  reported,  only  a  brief  sum¬ 
mary  will  be  given  here. 


The  equations  are  solved  in  dimensionless  parabolic  coordinates, 

4  and  ri,  which  are  related  to  the  dimensional  Cartesian  coordinates, 
x  and  y,  as  follows: 


x  +  iy 


£  +  it) 


where  xj  is  a  typical  distance  from  the  leading  edge,  and  Rx  ^  is  the 
Reynolds  number  based  on  the  mean  freestream  velocity  and  x^. 

The  time  is  made  dimensionless  with  the  mean  freestream  velocity 
and  with  x  j 


T  =  tU  /x  . 
oo  1 


(3) 


The  dimensionless  stream  function  and  vorticity  are  defined  in  terms 
of  the  dimensional  quantities  ^  and  £ . 


0  =  (2l/U  x  )1/2  $f  =  (2l/U  x,)1/2g 

OO  1  ^  oo  1  ° 


£  =  UfU  / (2fx . ) 

oo1  OO  1 


1/* 
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(4) 

(5) 
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where  V  is  the  kinematic  viscosity.  The  dependent  variable  f  is  the 
usual  Blasius  stream  function. 


The  equations  solved  are  the  parabolized  vorticity  equations  [8], 

*8,C '„/«  -8„(C/0{  +fB„rl/(2{2R  ) 

L  *'J« 

££  =  g  +  g  //ZR  \. 

^  xlj 


The  equation  set  is  solved  in  the  space 


1S?  sK/s)‘‘ 


0  ^  T)  <  oo 
0  <  T 

The  equations  are  first  order  in  time  and  are  integrated  from  a 
Blasius  initial  condition  at  r  =  0 


g  =  =  *f- 


5  s  B lasius 

The  four  boundary  conditions  in  the  p -dimension  are 

g  =  =  0;  P  =  0  (12) 

g  /£  =  f-p  =  1  +  fcostjt 

p  -*  oo  (13) 

g/£  =  f  =  P  -/9  +  £  (n  -y )  cos  ut 

where  yS  is  a  constant  characteristic  of  the  displacement  thickness 
and  y  is  either  yS/2  or  zero  in  order  to  be  consistent  with  the  bound¬ 
ary  conditions  at  Rx^,  discussed  subsequently. 

Posing  the  problem  on  a  finite  Reynolds  number  range  [re  Eq.  (8)] 
has  the  advantage  that  the  difficult  to  treat  vortical  singularity  at 
the  leading  edge  of  the  plate  is  avoided.  On  the  other  hand,  it  intro¬ 
duces  a  requirement  for  boundary  conditions  to  be  imposed  at 
The  boundary  conditions  are  obtained  from  the  unsteady  boundary- 
layer  solutions  [5-7].  The  small  s  =  cjx/Uqq  solution  for  the  x-com- 
ponent  of  the  velocity  is 

u  =  fg  +  £  cos(jt  [fg  +  pfg/2  +  O(s)]  (14) 

where  fg  is  the  Blasius  solution.  The  solution  for  large  s  is 

u  =  f'  +  £  £c ostut  -  exp(- ky)cos  (tut  -  ky)  +  0(s’^^)]  (15) 


where 


k  =  (u»/ 2.V) 


The  mixture  of  independent  variables  in  Eq.  (15)  is  used  to  empha¬ 
size  the  physics  of  that  equation.  Equation  (14)  shows  that  near  the 
leading  edge  the  unsteady  velocity  profile  is  in  phase  with  the  free 
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stream,  and  the  thickness  grows  at  the  same  rate  as  the  Blasius 
boundary-layer  thickness.  For  s  >  1,  Eq.  (15)  shows  that  the  un¬ 
steady  shear  layer  no  longer  grows  with  x,  and  that  there  is  a  phase 
shift  across  the  Stokes-type  layer.  For  large  s,  in  the  words  of 
Ackerberg  and  Phillips  |7J,  there  is  a  "double  boundary  layer"  with 
the  Stokes  layer  being  thinner  than,  and  embedded  in,  the  Blasius 
boundary  layer. 


I'he  upstream  boundary  conditions  obtained  from  Eqs.  (14)  and  (15) 
a  re 


K  r  g> 
s 


‘B 


(' 


B 


B 


(17) 


Equations  (17)  tollow  from  Eq.  (14);  the  conditions  obtained  from 
Eq.  (15)  are 


rjcoscot  +  |  exp( -Knlfcos  (wt  -  Kn)  +  sin(cot  -  Kp)J 
-  costot  -  sintjt  |  /lK 


«  r  *B  ** 


-  f  q  +  fn[cosiot  -  exp( -Kri  )cos  (cjt  -  Kp)] 


where  K  is  a  dimensionless  k 


K  =  kxl  Rx 


1/2 


(18) 


(19) 


The  parabolized  system  of  equations  is  third  order  in  £  and  requires 
a  downstream  boundary  condition.  The  present  numerical  results 
confirm  previous  experience  [8]  in  which  the  condition 

Sttt  =  0  (20) 

s  ss 

produced  a  small  disturbance  at  the  downstream  boundary. 

The  equations  are  solved  using  a  spectral  method  in  space;  a  Cheby- 
shev  polynomial  expansion  is  used  in  each  dimension.  The  solution 
is  updated  in  time  with  an  explicit  finite  difference  method. 

III.  Continuous  Wave -Wave  Interaction 

Two  possible  mechanisms  for  feeding  energy  into  a  Tollmien- 
Schlichting  wave  from  a  sound  wave  are  investigated  herein.  One 
possibility  is  the  existence  of  a  continuous  flow  of  energy  into  the 
Tollmien-Schlichting  wave  along  the  whole  boundary  layer.  The 
other  possibility  is  that  the  sound  wave  interacts  with  the  boundary 
layer  only  in  localized  regions,  and  then  the  Tollmien-Schlichting 
wave  propagates  independently  of  the  sound  wave.  In  this  section, 
the  former  possibility  is  investigated. 

I  he  study  of  the  wave-wave  interaction  is  complicated  by  the  fact 
that  the  computer  code  generates  solutions  for  a  finite  region  on  the 
plate.  Thus,  any  solution  obtained  is  dependent  upon  the  assumed 
upstream  boundary  conditions.  Furthermore,  unless  the  upstream 
boundary  conditions  are  specifically  chosen  to  be  orthogonal  to  the 
primary  eigensolution  of  the  Or r-Sommerfe Id  equation,  a  finite 
Tollmien-Schlichting  wave  will  be  generated  by  the  conditions.  In 
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this  section,  the  objective  is  to  investigate  the  possibility  of  a  con¬ 
tinuous  interaction.  This  continuous  interaction,  if  it  exists,  will 
be  most  easily  found  if  the  Tollmien-Schlichting  wave  generated  by 
the  upstream  boundary  condition  is  minimized.  As  will  be  shown 
subsequently,  Eqs.  (18)  generate  a  very  small  Tollmien-Schlichting 
wave.  The  effect  of  the  boundary  condition  on  the  solution  is  also 
clarified  by  moving  the  position  of  the  boundary  and  by  running  com¬ 
putations  in  the  Reynolds  number  range  in  which  the  boundary¬ 
generated  wave  decays. 

Before  the  solutions  are  examined,  some  trigonometric  identities 
will  be  presented  to  aid  in  their  interpretation.  If  at  some  station  in 
the  boundary  layer  the  sound  wave  generates  a  time  periodic  dis¬ 
turbance  of  amplitude  A  and  the  Tollmien-Schlichting  wave  generates 
a  disturbance  of  amplitude  B  and  real  wave  number  Of,  then  the  total 
disturbance  (at  frequency  u)  is 

u  =  A  cosut  +  B  coslax  -  (Jt)  (21) 

where  A  and  B  may  be  functions  of  y  and  slowly  varying  functions  of 
x,  and  the  origin  of  x  is  assumed  to  be  such  that  the  phase  angle  is 
zero.  The  Fourier  amplitude  of  the  disturbance  at  frequency  u  is 

2  2  1/2 

C.  =  ( A ~  +  B  +  2A  BcosOtx )  '  (22) 

When  A  »  B,  Eq.  (22)  may  be  expanded  to  give 

C  =  A  +  B  cosax  +  0(BP.  (23) 

Equation  (2  1)  assumes  that  A  and  B  are  independent  or  weakly  depen¬ 
dent  parameters.  Equation  (23)  shows  that,  in  this  case,  the  Fourier 
amplitude  of  the  total  disturbance  will  oscillate  spatially  with  a  wave 
length  equal  to  the  Tollmien -Sch hchting  wave  length,  that  the  magni¬ 
tude  of  the  spatial  mean  is  associated  with  the  sound  wave,  and  that 
the  magnitude  of  the  envelope  about  the  mean  is  associated  with  the 
Tollmien-Schlichting  wave.  The  experimental  data  of  Shapiro  [l]  are 
consistent  with  this  picture  in  that  he  found  the  Fourier  amplitude  to 
oscillate  spatially  with  the  Tollmien-Schlichting  wave  length. 

The  first  few  figures  show  the  results  of  a  numerical  calculation  in 
the  range  1.2  x  1  CP  -  Rx  —  3.  8  x  1(P.  The  dimensionless  frequency 
is  that  of  Shapiro,  uV/  =  5.  6  x  It)--.  The  large  s  boundary  layer 
equations  are  used  to  generate  the  upstream  boundary  conditions 
|Eqs.  (18)]  imposed  at  Rx  =  1.2  x  1  CP  or  s  =  6. 72.  The  nume  rical 
solutions  are  obtained  by  running  the  calculations  until  the  solution  is 
time  periodic;  this  solution  can  then  be  studied  by  Fourier  transform¬ 
ing  in  time . 

In  Fig.  1,  the  Fourier  amplitude  of  u  at  frequency  CJ  is  plotted  versus 
Reynolds  number  for  y  =  constant  ( Ry  =  yU^/v  ~  196).  This  value  of 
y  corresponds  to  tl  =  0.4  at  Rx  =  1.2  x  10-.  For  the  computation 
shown,  the  value  of  €  was  0,001  [re  Eq.  (1)].  (Numerical  evidence 
obtained  by  repeating  selected  computations  with  C  halved  verified 
that  solutions  with  this  value  of  e  are  linear  in  £.)  The  solution  has 
the  features  of  the  sum  of  a  Stokes  layer  solution  and  a  Tollmien- 
Schlichting  wave  solution.  The  mean  curve  in  Fig.  1  is  very  nearly 


ant.  corresponding  to  the  fact  that  the  high  Reynolds  number 
ady  boundary  layer  solution  (Stokes  layer!  is  a  function  only 
The  solution  oscillates  about  the  mean  with  a  wave  length  equal 
■  Tollmien  -  Schlichting  wave  length  (shown  as  A  in  T  ig.  1). 

ler  to  confirm  the  presence  of  a  Tollmien- Schlichting  wave,  the 
ttude  ot  the  envelope  of  several  curves  similar  to  tig.  1  was 
net!  and  plotted  versus  n  m  Tig.  2  at  Rx  2.8  x  lt)\  The  points 
tamed  are  compared  with  a  scaled  solution  to  the  O  r  r  -  Somme  r  fe  Id 
ion  at  the  same  Reynolds  number,  and  the  agreement  is  excellent. 
*e  2  shows  that  a  Tollmien  -  Schlichting  wave  is  present,  anil  that 
nvelope  of  curves  such  as  T  ig.  1  characterize  the  magnitude  of 
a  ve . 


He  envelope  points  are  shown  for  values  of  n  greater  than  about 
because  at  the  larger  values  of  n  the  solution  curves  appear  to  be 
la  ted  with  two  spatial  wave  lengths  rather  than  just  one,  making 
xtraction  of  the  required  information  very  difficult.  Figure  3 
s  a  plot  of  such  a  solution  curve  at  a  location  twelve  times  further 
the  wall  than  the  solution  in  T  ig.  1.  The  value  of  y  in  Fig.  3 

•  sponds  to  a  value  of  ri  of  4.  8  at  Rx  1.  2  x  10*1  and,  for  com  pari 
.ith  t  ig.  2,  n  %  ' ,  1  1  at  =  2.  8  x  10  '.  The  spatial  modulation 

•  Tollmien  Schlichting  wave  length  is  still  in  evidence,  but  a 

vith  a  .  it  three  tin  es  as  long  is  also  present, 

•ms  likely  that  the  longer  wavelength  modulation  indicates  the 
•n  e  >t  listurbances  described  by  the  continuous  spectrum  of  the 
S  mm  e  t  ie  Id  equation  (.  ontinuous  spectrum  solutions  are  very 

1  m  the  boundary  layer  ompared  to  their  value  further  from  the 
I  hi'  is  onsistent  with  the  fact  that  the  higher  wave  length 
ilation  is  not  seen  >n  !  ig.  1.  Furthermore,  the  dimensionless 
speed  !  the  continuous  spectrum  modes  is  very  close  to  unity, 
iared  t.  a  lollr  ien  -  Schlichting  wave  speed,  in  the  present  case, 

.  Thus,  the  factor  of  three  difference  in  wave  speed  is 
stent  with  the  factor  >f  three  difference  in  modulation  wave  length 
g.  3. 
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The  main  interest  here  is  the  generation  of  Tollmien-Schlichting 
waves;  thus,  most  of  the  remainder  of  the  plots  will  be  at  a  y- loca¬ 
tion  very  near  the  wall  (Ry.  =  196).  In  all  cases,  this  was  sufficient 
to  suppress  the  continuous  spectrum  disturbances  which  were  evi¬ 
dent  further  out  in  the  boundary  layer.  Any  other  discrete  eigenmodes 
present  decay  very  rapidly. 

Figure  4  shows  the  variation  of  the  total  Fourier  amplitude  of  u  versus 
y  at  the  upstream  boundary,  and  in  the  interior  of  computation  domain. 
The  upstream  boundary  condition  has  been  imposed  at  Rx  =  1.  2  x  10^, 
so  the  left  curve  in  Fig.  4  is  just  the  Fourier  amplitude  of  Eq.  (15). 
This  figure  demonstrates  that  the  high  Reynolds  number  boundary  layer 
solution  is  an  excellent  approximation  to  the  solution  of  the  more  com¬ 
plete  equations.  Any  deviation  of  the  right-hand  curve  of  Fig.  4  from 
the  (left-hand)  bounda ry- layer  solution  because  of  the  presence  of  the 
Tollmien-Schlichting  or  continuous  spectrum  waves  is  imperceptible 
on  the  scale  of  Fig.  4. 


Fig.  3.  Fourier  amplitude  at 
constant  y 


y(Uoc/2i/x1)1/2 

Fig.  4.  Fourier  amplitude  of 
u  versus  y 


Figures  1  through  4  demonstrate  that  a  Tollmien-Schlichting  wave  is 
present.  The  main  interest  of  this  paper  is  the  source  of  this  wave. 
The  spatial  amplitude  variation  of  the  Tollmien-Schlichting  wave  in 
the  presence  of  the  sound  wave  is  compared  with  the  variation  of  a 
pure  Tollmien-Schlichting  wave  with  no  sound  wave  present.  If  the 
curves  have  the  same  behavior,  then  the  Tollmien-Schlichting  wave 
propagates  independently  of  the  sound  wave;  if  the  envelope  curve 
grows  relative  to  the  pure  wave  amplitude,  then  energy  is  fed  into  the 
Tollmien-Schlichting  wave  continuously  by  the  sound  wave.  The  cir¬ 
cles  in  Fig.  5  are  obtained  from  the  extrema  of  Fig.  1;  the  curve  is 
a  scaled  linear  Tollmien-Schlichting  wave  solution  with  no  sound  wave 
present.  With  the  exception  of  the  first  point,  they  are  in  good  agree¬ 
ment.  At  the  first  point,  other  decaying  eigenmodes  are  non-negli- 
gible.  It  is  concluded  that  the  Tollmien-Schlichting  wave  in  the  nu¬ 
merical  computation  is  generated  by  the  boundary  condition  at  the 
upstream  boundary  and  subsequently  propagates  independently  of  the 
sound  wave.  As  noted  previously,  Thomas  and  Lekoudis  [2]  have 
concluded  from  Shapiro's  data  that  the  two  waves  are  independent. 
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Figure*  (■  further  demons!  rates  that  the  Tollmien  Schliehting  w.n-e 
amplitude  in  the  present  eomputat  ions  is  eont  rolled  by  the  upstream 
boundary  eonditions.  1  lie  upstream  boundary  has  been  moved  from 
l .  x  10  in  V  i|t.  V  to  1.  I  \  10  '  in  Fig.  t>;  all  other  parameters  are 
unchanged.  The  phase  and  amplitude  are  signi  Meant  lv  atteeted  l»v  the 
change  in  location  ot  the  upstream  boundary,  rather  than  being  indo 
pendent  ot  the  boundary  condition  as  would  be  the  ease  if  interaction 
with  the  treestream  sound  wave  dominated. 


VO  s  Kv  10  S  Kv 


Fig.  •**.  i  omparison  of  envelope  Fig.  ti.  Fourier  amplitude  of 

behavior  with  a  pure  Tollmien  u  at  constant  y 

Sell  holding  wave 

1  he  main  result  of  this  section  is  that  a  small  amplitude  rollmion 
Schliehting  wave  may  be  added  tv*  the  Stokes  solution  and  that  these 
wa  ves  iir  solutions  do  not  interact,  in  the  high  Reynolds  number  range 
considered.  Hie  rollmion  Schliehting  waves  present  m  the  solutions 
discussed  were  found  tv'  he  generated  artificially  by  the  upstream 
bound  a  rv  condit  ion . 

H  .  Focalized  Wave  Wave  Interaction 

In  the  previous  section,  it  was  concluded  that  the  interaction  between 
a  sound  wave  and  a  I'ollmien  Schliehting  wave  does  not  occur  over  the 
entire  boundary  layer  but,  it  it  occurs,  the  interaction  is  confined  to  a 
localized  region.  1  he  purpose  ot  this  section  is  to  determine  where 
energy  is  ted  into  the  I'ollmien  Sell  It  cht  ing  wave  and  to  compute  the 
amplitude  ot  this  wave. 

Some  ot  the  major  results  ot  unsteady  boundary  layer  theory  will  he 
summarized  in  order  to  place  the  numerical  results  in  perspective. 

In  the  boundary  layer  rase,  there  are  two  distinct  solution  regions  on 
tin*  plate.  At  low  I  \evnolds  number,  the  unsteady  layer  grows  at  the 
same  rate  as  the  lllasius  boundary  layer;  at  high  Reynolds  number,  the 
unsteady  layer  is  controlled  by  the  frequency  and  is  thinner  than  the 
lllasius  layer.  There  i  s  little  doubt  that  the  more  complete  equations 
considered  herein  have  solutions  with  essentially  the  s.  me  behavior  as 
Fig.  *1  has  demonstrated.  The  authors  ot  Rets,  o  and  7  were  concerned 
with  explaining  how  the  parabolic  boundary  layer  equations  could  al 
ways  have  the  same  large  Reynolds  number  solution,  independent  ot 


.my  di  sturbance  s  which  might  exist  upstream.  They  conclude  that  the 
downstream  solution  is  approached  via  damped  eigenfunction* ;  dis¬ 
turbances  inconsistent  with  the  downstream  solution  generate  asymp¬ 
totic  eigensolut  ions  which  decay  exponentially.  Numerical  solutions 
to  the  parabolized  vorticity  equations  will  be  shown  to  be  similar. 

Even  with  large  disturbances  (introduced  via  the  upstream  boundary 
conditions),  the  solution  at  large  Reynolds  number  approaches  the 
expected  downstream  solution  via  eigenfunctions.  There  is  one  major 
difference,  however;  the  asymptotic  eigenfunction  (not  present  in  the 
boundary  layer  approximation)  associated  with  the  Tollmien -Schlichting 
wave  is  not  necessarily  damped.  The  following  picture  will  be  de¬ 
veloped.  Under  normal  flow  conditions  (no  localized  disturbance),  the 
magnitude  of  the  various  eigenfunctions  including  the  Tollmlen- 
Schlichting  wave  amplitude  in  the  case  of  the  more  complete  equations 
is  set  by  the  end  of  the  region  in  which  the  unsteady  bounds  ry  layer 
changes  from  a  growing  to  a  constant- thickness  layer.  If  disturbances 
are  introduced  at  some  point  downstream  of  this  region,  then  the 
magnitude  of  the  eigenfunctions  is  set  immediately  downstream  of  the 
disturbance. 

All  the  numerical  results  reported  in  this  section  use  Eq.  (17)  as  the 
upstream  boundary  condition  on  the  computation.  Physically,  the  use 
of  this  boundary  condition  is  only  justified  for  s  small  compared  to 
unity.  Computations  were  actually  performed  in  which  the  upstream 
value  of  s  was  varied  from  It'  to  0.  Using  Eq.  (17)  at  high 

Reynolds  number  is  a  way  of  inputting  a  large  disturbance  at  high  Rev 
tiolds  number.  Results  from  these  computations  will  be  used  to  demon 
strate  that  the  solution  does  relax  to  the  expected  high  Reynolds  mini 
her  behavior. 

Numerical  instabilities  were  encountered  when  the  first  attempts  were 
made  at  running  computations  with  small  upstream  Reynolds  number. 

It  was  ultimately  found  that  with  ''  Chebyshev  modes  in  the  stream 
wise  direction,  the  requirement  for  numerical  stability  is 

R  /R  <  F  '>  (24) 

xi 

this  is  ,i  severe  restriction,  because  the  physical  length  ot  the  com 
putation  varies  linearly  with  the  upstream  Reynolds  number  Rx  j .  and 
the  interest  is  in  small  values  of  Rx..  At  larger  values  of  Rx  ,  the 
length  of  the  computation  is  restricted  to  about  eight  I'ollmien  1 
Schlichting  wavelengths  with  the  '  t  modes.  Uecause  of  the  restriction, 
no  computations  were  made  with  a  value  ot  s  at  the  upstream  boundary 
smallrr  t  ha  n  0 .  t  '>2 . 

Consider  first  a  calculation  which  is  identical  tv'  Uig.  except  that 
the  upstream  boundary  condition  has  been  changed  from  Eq.  (Is)  tv' 

Eq.  (14).  Vhe  boundary  condition  is  physically  correct  for  small  s. 
so  it  is  here  interpreted  as  a  solution  with  a  large  periodic  disturbance 
input  at  s  *  ('.  It',  Rx  =  1.1x10'.  A  plot  of  the  Fourier  amplitude 
of  ii  at  constant  v  analogous  to  Fig.  <-  is  shown  in  Fig.  7.  The  two 
curves  are  qualitatively  very  similar;  after  a  finite  region  at  the  left- 
hand  etui  in  which  continuous  ami  higher  discrete  modes  decay,  both 
are  modulated  with  the  same  I’ollmien  Schlichting  wavelength  a  ml 
have  the  same  envelope  decay  rate  and  mean  value.  The  major  dif 
terence  is  that  the  I'ollmien  Schlichting  wave  (the  envelope)  amplitude 


•  i 


Fourier  amplitude  of  u/« 


1CT5  Rv 


Fig.  7.  Fourier  amplitude 
of  u  at  constant  y 


is  about  a  factor  of  16  larger  in  Fig.  7  than  in  Fig.  6.  It  is  also 
interesting  to  Look  at  the  profile  of  the  Fourier  amplitude  of  u  versus 
y  and  see  how  it  changes  as  a  function  of  Reynolds  number.  Such 
curves  are  shown  in  Figs.  8  and  9.  The  curve  for  Reynolds  number 
1.  1  x  10^  in  Fig.  8  is  the  upstream  boundary  condition  or,  alterna¬ 
tively,  the  unsteady  disturbance.  The  characteristic  thickness  of  this 
disturbance  is  much  greater  than  the  normal  Stokes  layer  which  would 
exist  at  that  frequency  and  Reynolds  number.  Figures  8  and  9  show 
that  the  boundary  layer  has  adjusted  to  a  solution  similar  to  those 
shown  in  Fig.  4  by  Reynolds  number  2  x  10^,  Downstream  of 
Rx  =  2  x  1(P  the  characteristic  shape  shown  in  Fig.  4  persists,  sub¬ 
ject  to  some  small  modification  near  the  peak  of  the  Tollmien- 
Schlichting  wave  (ri  ~  11.  Figures  6-8  show  that  the  present  solutions 
have  the  same  behavior  as  the  boundary  layer  solutions  - -they  approach 
the  same  high  Reynolds  number  solution  via  their  respective  eigen¬ 
functions. 

It  was  not  possible  to  carry  out  numerical  computations  with  the  up¬ 
stream  boundary  at  an  arbitrarily  small  Reynolds  number.  Thus, 
computations  have  been  carried  out  with  several  upstream  boundary 
locations,  and  the  results  are  extrapolated  to  smaller  values. 

In  Fig.  10,  the  amplitude  of  several  Tollmien -Schlichting  waves 
generated  by  imposing  the  boundary  condition  (14)  at  different  loca¬ 
tions  is  plotted  versus  Reynolds  number.  These  curves  are  obtained 
from  the  respective  envelope  curves,  the  uppermost  curve  being 
obtained  from  Fig.  7.  There  is  no  numerical  evidence  that  the  ampli¬ 
tude  of  the  Tollmien-Schlichting  waves  is  affected  by  the  presence  of 
the  sound  wave  subsequent  to  their  generation  at  the  upstream  bound¬ 
ary.  This  may  be  verified  by  noting  that  to  within  numerical  accuracy 
the  curves  in  Fig.  10  are  each  a  constant  multiple  of  any  other  curve. 
This  is  the  expected  behavior  of  a  linear  Tollmien-Schlichting  wave; 
the  rate  of  decay-growth  is  independent  of  amplitude.  In  those  cases 
in  which  the  computation  did  not  extend  to  the  neutral  (minimum)  point, 
Rx  «  3  x  105,  the  curves  have  been  extrapolated  using  the  appropriate 
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Fig.  8.  Variation  of  Fourier  Fig,  4.  Variation  of  Fourier 

amplitude  across  boundary  amplitude  across  boundary 

layer  layer 

constant  and  are  shown  on  the  figure  as  dashed  lines.  Each  of  th»> 
curves  in  Fig.  10  was  obtained  from  a  single  numerical  computation 
except  for  the  lower  curve  with  Rx  ^  -  10"*.  This  curve  was  obtained 
with  three  overlapping  computations  with  Reynolds  number  ranges 
104  <  Rx  <  l.  9  x  10-*,  2Q-187  <Hs<  1.  IS  x  105,  and  82h2t>  <  R  < 

'.2  x  10s.  Restarting  the  calculation  introduces  a  substantial 
amount  of  noise  into  the  computation.  Because  of  this  noise  problem 
and  the  limit  in  F.q.  (241,  no  computations  were  attempted  with  Rx 
less  than  7  x  10.  Increasing  the  number  of  modes  in  the  £  -direct 
tion  would  allow  Rx  to  be  reduced  further  but  was  rejected  because 
of  computer  cost  considerations. 

The  numerical  results  presented  are  consistent  with  the  boundary 
layer  model  studied  in  Refs.  6  and  7.  In  the  boundary  layer  case, 
the  large  Reynolds  number  solution  is  approached  via  asymptotic, 
decaying  eigenfunctions.  r he  amplitude  of  these  eigenfunctions  is 
determined  uniq  lely  bv  the  upstream  boundary  condition,  and  the 
sound  wave  does  not  feed  further  energy  into  these  eigenmodes  at 
downstream  locations.  In  the  present  case,  the  eigenmodes  of  interest 
are  those  associated  with  the  Or r -Sommer fe Id  equation  and  are  dif¬ 
ferent  from  the  modes  of  the  boundary  layer  equations.  The  numerics 
have  been  concentrated  on  the  Tollmien -Schlichting  mode.  Analogous 
to  the  bounda  ry  -  laye  r  case,  the  amplitude  of  this  mode  is  set  by  the 
upstream  boundary  condition,  and  no  further  energy  is  found  to  be  fed 
into  it  at  downstream  stations.  In  the  limiting  case  wherein  the 
boundary  condition  is  imposed  arbitrarily  close  to  the  leading  edge, 
the  eigenmodes  in  both  cases  are  still  uniquely  determined.  How¬ 
ever.  these  eigenmodes  are  asymptotic  solutions  valid  only  for  s  >  l. 

(It  is  shown  in  the  next  section  that  the  solutions  to  the  Orr-Sommer- 
fe  Id  equation  only  make  physical  sense  for  s  >  1.1  In  this  context,  the 
region  OS  8  <  (1(11  is  the  region  in  which  the  boundary  layer  and  the 
sound  wave  interact.  The  output  of  this  region  is  a  set  of  asymptotic 
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Fig.  10.  Fourier  envelope  of  u  Fig.  11.  Amplitude  at  neutral 
versus  Reynolds  number  point  versus  R 


linear  solutions  which  are  physically  and  mathematically  independent. 
Among  these  independent  solutions  are  the  Tollmien-Schlichting  wave 
and  the  Stokes  layer  solution. 

The  neutral  (minimum)  amplitude,  obtained  from  Fig.  10  and  similar 
curves,  is  plotted  in  Fig,  11  as  a  function  of  the  upstream  boundary 
location.  l'his  includes  three  points  obtained  from  computations  with 
s  <  1  at  RXj.  The  neutral  value  for  Rx  j  0  is  the  amplitude  of  a 
1  ollmien  -  Schlichting  wave  generated  by  a  sound  wave  with  dimension- 
tc-s  frequency  51  x  F)  "  at  R  «  ll>t,.  By  extrapolation  of  Fig.  11, 
this  neutral  amplitude  is  estimated  to  be  5.5  *  1 .  5  x  10~Y  The 
neutral  point  amplitude  of  the  Tollmien  Schlichting  wave  at  R  =  1% 
f  Fig.  2 )  is  about  half  the  peak  amplitude.  Thus,  the  peak  amplitude 
of  the  sound  generated  lollmien-Schlichting  wave  is  1.0  ±  0.3  x  10~4 
times  the  sound  wave  amplitude. 

1  he  frequency  of  the  sound  wave  and  the  lollmien-Schlichting  wave  in 
the  computations  described  up  to  this  point  is  =  5t>  x  10'“. 

l’his  value  was  chosen  to  coincide  with  Shapiro's  |  1^  experimental 
frequency.  There  is  substantial  disagreement  between  this  experi¬ 
mental  work  and  the  present  numerical  computations.  Based  upon 
Shapiro's  Figure  28,  the  Tollmien-Schlichting  wave  amplitude  at  the 
neutral  point  is  about  0.4  times  the  sound  wave  amplitude,  while  the 
present  computations  suggest  a  value  of  1()'4.  There  are  two  possible 
reasons  tor  this  discrepancy.  The  data  were  obtained  on  an  elliptical¬ 
nosed  model  with  a  significant  pressure  gradient  up  to  Reynolds  num¬ 
bers  greater  than  10s.  Thus,  the  data  and  the  computations  are  not 
strictly  comparable.  Another  possibility  is  that  the  sound  wave  in  the 


experiments  was  not  exactly  a  plane  wave  and  thus  effectively  gene¬ 
rated  local  disturbances.  The  numerical  results  of  this  section  sug¬ 
gest  that  a  localized  disturbance  can  generate  a  larger  Tollmien- 
Schlichting  wave  than  a  plane  wave  disturbance. 

A  brief  investigation  has  been  made  with  a  smaller  frequency  of  the 
disturbing  sound  wave.  Figure  12  is  a  plot  identical  to  Fig.  10, 
except  that  the  dimensionless  frequency  has  been  halved  to  28  x  10'^. 
This  frequency  is  approximately  the  frequency  at  which  the  total 
growth  of  the  Tollmien-Schlichting  wave  in  the  unstable  region  is  e  ', 
as  shown  for  example  by  Jaffe,  Okamura  and  Smith  [  1 0 ) .  Although 
only  two  computations  have  been  made  at  this  frequency,  the  be¬ 
havior  is  similar  to  that  at  the  higher  frequency.  The  neutral  point 
amplitudes  from  Fig.  12  are  shown  in  Fig.  13.  A  linear  extrapola¬ 
tion  of  these  results  to  zero  upstream  Reynolds  number  gives  a  wave 
amplitude  of  3  x  1 0 _ ^  times  the  sound  wave  amplitude.  Note,  how¬ 
ever,  that  the  results  in  Fig.  11  suggest  that  the  linear  extrapolation 
may  be  a  factor  of  two  or  three  too  high.  Again,  taking  into  account 
the  fact  that  the  plotted  results  are  not  at  the  peak  of  the  Tollmien- 
Schlichting  wave,  the  linearly  extrapolated  result  gives  a  peak  wave 
amplitude  at  the  first  neutral  point  of  7  x  10*5  times  the  sound  wave 
amplitude.  Thus,  it  is  concluded  that  halving  the  sound  wave  frequency 
decreases  the  peak  Tollmien-Schlichting  wave  amplitude  by  30  to  70 
percent. 


On  a  sharp,  flat  plate,  a  plane  sound  wave  interacts  with  the  boundary 
layer  very  near  the  leading  edge  and  generates  a  Tollmien-Schichting 
wave  which  can  first  be  identified  as  such  downstream  of  the  s  !s  0(1) 
region.  This  wave  propagates  unaffected  by  the  sound  wave  after  it 
appears.  Between  s  =  0(1)  and  the  neutral  point,  the  wave  decays 
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Fig.  12.  Fourier  amplitude  of 
u  versus  Reynolds  number 


Fig.  13.  Amplitude  at  neutral 
point  versus  R^ 
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substantially,  to  about  10  4  times  the  sound  wave  amplitude  at  the 
neutral  point.  Ihis  is  small,  but  it  should  he  noted  that  a  total 
Tollmien -Sch  lie  h  ting  wave  growth  of  e9  =  104  downstream  of  the 
neutral  point  is  frequently  used  to  empirically  correlate  transition. 
The  results  of  this  section  also  suggest  that  localized  disturbances 
may  generate  larger  amplitude  Tollmien-Schlichting  waves  than  the 
plane  wave  disturbance. 

\  .  Solutions  to  the  Orr -Somme  rfeld  Equation 

In  this  section,  the  behavior  of  solutions  to  the  Or r -Somme r fe Id 
equation  in  the  regions  s  <  0(1)  and  s  >  0(1)  are  studied.  These 
solutions  provide  further  confirmation  of  the  difference  in  the  phy 
sics  ot  unsteady  boundary  layers  in  the  two  regions. 

The  Or  r-Somme rfeld  equation  (i.e.,  see  Schlichting  [11])  is  obtained 
by  linearizing  the  Na  vie  r -Stokes  equations  and  then  representing  the 
disturbance  in  the  form  0(y)  exp  i(ctx-<Jt).  The  spatial  eigenfunc¬ 
tions  in  which  u  is  real  andot  is  complex  are  of  interest  here.  It 
should  be  emphasized  that  the  variation  of  Of  with  x  or  R  is  implicitly 
assumed  to  be  negligible,  since  the  derivative  with  respect  to  x  is 
taken  as  iCt/>exp  i(OI\  ut). 


I  he  primary  spatial  eigenvalue  is  plotted  in  Fig.  14  as  a  function  of 
Reynolds  number.  The  emphasis  is  on  low  Reynolds  number,  and 
the  well-known  unstable  region  is  not  computed.  (The  numerical 
method  used  is  described  in  the  Appendix  of  Ref.  9.)  Above  a  Reynolds 
number  of  about  6.5  x  104,  the  real  part  of  a  has  a  slight  negative 
slope  corresponding  to  the  small  increase  in  Tollmien-Schlichting 
wave  length  with  Reynolds  number.  Starting  from  the  continuous  spec- 
trum  at  R  ~  62.9  Reap/ 1’  grows,  rapidly  at  first,  from  a  value  of 
about  0.  _  *•>  x  10-  .  (The  relationship  between  the  discrete  and  con¬ 
tinuous  spectrum  is  discussed  in  Refs.  9  and  12.)  The  imaginary  part 
Ot  at'  I  OH  is  about  2)2  x  10*  at  the  continuous  spectrum  and  decreases 
by  two  orders  of  magnitude  with  the  higher  Reynolds  number  behavior 
as  shown  in  I'ig.  14. 


I  he  minimum  wavelength  associated  with  the  maximum  of  Red  in 
1  ‘K\  *4  lf  KX  =  3-  *  *  1°  :  thus.  6e low  a  Reynolds  number  of  about 

'  *  -  i  !t-Cha.nR!u  l^a  °^'Ur  ln  °ne  wave  l^gth.  The  assumption 

°X  leading  to  the  Or  r-Somme  rfeld  equation  is  violated  and  a 

different  physical  model  is  required  in  this  region.  Thus,  the  physical 
entity  which  is  approximately  modeled  by  the  Or  r-Somme  rfeld  equa¬ 
tion  and  known  as  a  rollmien-Schlichting  wave  can  only  exist  above 
some  finite  Reynolds  number;  in  the  present  case  this  is  about  3  x  104. 


The  solutions  of  the  Orr -Somme  rfe  Id  equation  are  consistent  with  the 
numerical  solutions  of  the  more  complete  equations  in  that  two  distinct 
regions  are  found.  For  s  >  0(1),  where  the  Tollmien-Schlichting  and 
sound  waves  are  independent,  the  Or r-Sommerfeld  model  is  satisfac¬ 
tory.  For  s  <  0(1),  energy  is  fed  into  the  Tollmien-Schlichting  wave 
which  appears  downstream,  and  the  Ur  r  -  Somme  rfe  Id  equation  breaks 
down  as  a  valid  physical  model. 


1 

I 


Fig.  14.  Spatial  eigenvalues  of  the  Orr-Sommerfeld  equation 
VI.  Concluding  Hema  rk s 

The  major  results  to  emerge  from  this  study  are  as  follows.  For 
nearly  incompressible  flow  (sound  wave  length  large  compared  to 
'T'ollmien -Schlichting  wave  length),  plane  sound  waves  feed  energy  into 
To  llmien  -  Sch  lichting  waves  only  very  near  the  leading  edge;  over 
most  of  the  boundary  layer  the  two  waves  are  totally  independent.  Pre¬ 
dictions  have  been  made  of  the  amplitude  of  a  Tollmien-Schlichting 
wave  generated  by  tree-stream  disturbances;  the  physical  under¬ 
standing  achieved  in  this  work  will  aid  future  studies  of  other  types 
of  free -stream  distvirbances.  The  response  of  the  boundary  layer 
to  localized  unsteady  disturbances  has  been  extensively  studied. 

These  results  show  that  much  larger  Tollmien-Schlichting  waves  can 
be  generated  by  local  disturbances  and,  thus,  should  be  more  impor¬ 
tant  in  the  ultimate  transition  process.  Solutions  of  the  Orr-Sommer- 
teld  equation  are  found  to  change  character  near  s  =  0(1);  the  solutions 
for  s  <  0(1)  suggest  that  Tollmien-Schlichting  waves  do  not  occur  in 
this  small  region  near  the  leading  edge  of  the  plate. 
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